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The suppressed decay chain B~ — > DK~ , D — > K + ir~ , where D indicates & D or D state, pro- 
vides important information on the CP- violating angle fa. We measure the ratio IZdk of the decay 
rates to the favored mode B~ ->• DK~ , D -> K~tv + toheH D K = [1.63i^(stat)i£?£(syst)]xl(T 2 , 
which indicates the first evidence of the signal with a significance of 4.1o. We also measure the asym- 
metry Adk between the charge-conjugate decays to be Adk = — 0.39lg'2g(stat)lQ {^(syst). The 
results are based on the full 772 x 10 6 BB pair data sample collected at the T(4S) resonance with 
the Belle detector. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



Determinations of the parameters of the standard 
model arc fundamentally important; any significant dis- 
crepancy between the expected and measured values 
would be a signature of new physics. The Cabibbo- 
Kobayashi-Maskawa matrix 0, Q consists of weak inter- 
action parameters for the quark sector, one of which is 
the CP-violating angle fa = arg (~V ud V ub * /V cd V cb *) 0. 
Several methods proposed for measuring fa exploit in- 
terference in the decay B~ — > DK~ (D = D° or D° 
where the two D states decay to a common final state 
7|. One of the methods utilizes the decay B~ — > DK~ , 
D — > K + tt~ [||. The magnitudes of interfering ampli- 
tudes are comparable and hence can enhance the effects 
of CP violation. Previous studies of this decay mode 
have not found a significant signal yield 0, Q . 

In this analysis, we measure the ratio IZdk of the 
aforementioned suppressed decay to the favored decay, 
B~ — > DK~, D — > K~tt + , and the CP asymmetry Adk 
defined as 



TZ-dk 
■Adk 



B{[K + n-} D K-) + B{[K-tt + ]d K+) 
B([K-tt+] d K-) + B{[K+ir-] D K+) 
B{[K+ir-] D K-) - B{[K-tt+] d K+) 
B([K+tt-] d K-) + B{[K-ir+] D K+) 



, (1) 
, (2) 



where {f]o indicates the final state / originating from 
a D° or D° meson. The same selection criteria and fit- 
ting functions are used for the suppressed decays and the 



favored decays wherever possible in order to cancel sys- 
tematic uncertainties. The observables are related to fa 
as follows: 

IZdk = r% + r 2 D + 2r B r D cos (S B + S D ) cos fa, (3) 
Adk = 2r B rrj sin (5 B +S D ) sin fa/lZ dk, (4) 



where r B = \A(B~ -> D°K~)/A(B- -> D K~)\, 
r D = \A(D° -> K+ir-)/A(D° -> K+Tr~)\, and S B (S D ) 
is the strong phase difference between the two B (D) de- 
cay amplitudes appearing in the ratios. By combining 



other experimental inputs l3"12|, the value of 63 can 
be extracted in a model-independent manner Q . The 
decay B~ —> Dtt~, D — > K + it~ is also analyzed as a 
reference. For this final state the decay rate is relatively 
large whereas the effect of fa is small. 

The results are based on the full 772 x 10 6 BB pair data 
sample collected at the T(45) resonance with the Belle 
detector located at the KEKB asymmetric-energy e + e~ 
collider I13B . The Belle detector is described in detail 
elsewhere The primary detector components used 

in this analysis are a tracking system consisting of a sili- 
con vertex detector and a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), and a barrel-like arrangement of time-of- flight 
scintillation counters (TOF). 

Neutral D meson candidates are reconstructed from 
pairs of oppositely charged tracks. For each track, we 



apply a particle identification requirement based on in- 
formation from the ACC and TOF, and specific ioniza- 
tion measurements from the CDC. The efficiency to iden- 
tify a kaon or a pion is 85-95%, while the probability of 
misidentifying a pion (kaon) as a kaon (pion) is 10-20%. 
The invariant mass of the Kir pair is required to satisfy 
1.850 GeV/c 2 < M(Kir) < 1.880 GeV/c 2 , which corre- 
sponds to approxim ately ±3<r around the world-average 
value of the D mass [l5j , where a denotes the resolution 
in M(Ktt). To improve the momentum determination, 
tracks from the D candidate are refitted with their in- 
variant mass constrained to the nominal D mass. 

B meson candidates are reconstructed by combining a 
D candidate with a charged hadron candidate. The sig- 
nal is identified using the beam-energy-constrained mass 
(Mbc) and the energy difference (AE) defined, in the 
e+e - center-of-mass frame, as Mb c = y/ E^ cam — \pb\' 2 
and AE = Eb — £bcam, where -Ebeam is the beam en- 
ergy and ps and Eb are the momentum and energy, 
respectively, of the B meson candidate. We require 
5.271 GeV/c 2 < M bc < 5.287 GeV/c 2 , which corre- 
sponds to ±3cr around the B mass value [THJ] with a 
denoting the resolution in M bc . 

The dominant backgrounds arise from the continuum 
processes e + e~ — » qq (q = u,d,s,c). In order to re- 
move D*^ — > Dir^ decays produced in such a process, 
we employ the variable AM defined as the mass differ- 
ence between the D*^ and D candidates, where the D*^ 
candidate is reconstructed from the D candidate used in 
the B reconstruction and a tt^ candidate not used in the 
B reconstruction. We require AM > 0.15 GeV/c 2 , which 
removes 28% of the cc background as well as some of the 
BB background. The loss of signal efficiency is 1.4%. 

The qq background is further discriminated with a neu- 
ral network technique based on the NeuroBayes pack- 
age [l6j]. The inputs to the network are: 1) a Fisher 
discriminant 17| formed from modified Fox- Wolfram mo- 
ments [HI , 2) the vertex separation between the B candi- 
date and the remaining tracks, 3) the cosine of the decay 
angle of D — > K + ir~ , where the decay angle is defined as 
the angle between the K + candidate and the B~ candi- 
date in the rest frame of the D, 4) the cosine of the angle 
between the B candidate and the beam axis in the e + e~ 
center-of-mass frame, 5) the expected B flavor dilution 
factor that ranges from zero for no flavor information 
to unity for unambiguous flavor assignment (l9j . 6) the 
cosine of the angle between the thrust axis of the B can- 
didate and that of the rest of the event, where the thrust 
axis is oriented in such a way that the sum of momen- 
tum projections is maximized, and four other variables 
that exploit the kinematics of the events. The neural 
network is trained with Monte Carlo (MC) events. A re- 
quirement is applied on the network output {NB) that 
preserves 96% of the signal while rejecting 74% of the 
background. 

There are a few background modes that can peak in 



the signal window ('peaking background'). The decay 
B~ —> [K + K ~]dtt~ may contribute to the background 
for B~ — > [K + it~]dK~ if the D candidate is misrecon- 
structed. To reject this background, we veto candidates 
satisfying 1.840 GeV/c 2 < M (KK) < 1.890 GeV/c 2 . 
The favored decay B~ — > [K~ir + ]r>h~~ (h = K or it) can 
also produce a peaking background for the suppressed 
decay if both the kaon and the pion from the D decay 
are misidentified and the particle assignments are inter- 
changed. We thus veto candidates for which the invariant 
mass of the Kit pair is inside a (1.865 ± 0.020) GeV/c 2 
window when the mass assignments are exchanged. 

The signal yield is extracted from the two-dimensional 
distribution of AE and NB using an extended unbinned 
maximum likelihood fit. The fit is simultaneously applied 
to DK~ , DK + , Dir~ and Dn + . The components of the 
fit are divided into signal, 'feed-across', BB background, 
and qq background, as described in detail below. For 
each component, the correlation between AE and NB 
is found to be small. We thus obtain the probability 
density function (PDF) by taking a product of individual 
PDFs for AE and NB. For NB, we use one-dimensional 
histogram PDFs for all components. 

For the AE signal PDF, we use a sum of two Gaussians 
whose parameters are fixed from the data for the favored 
modes. For NB, we obtain the PDF by applying \AE\ < 
0.01 GeV to the same samples. 

The Dir (DK) feed-across is the background from 
misidentified Dir {DK) final states that peaks in the fit 
to the DK (Dir) sample. The shift due to the incorrect 
mass assignment makes the AE distribution asymmet- 
ric, and thus we use a sum of two asymmetric Gaus- 
sians, for which the left and right sides have different 
widths. The corresponding parameters for the Dir feed- 
across are obtained from the favored mode in data, while 
those for the DK feed-across are obtained from MC sam- 
ple because of low statistics in the data. The PDFs for 
NB are obtained from the same reference samples used 
for AE calibration after applying the additional require- 
ment \AE - 0.05 GcV| < 0.01 GeV to the data sample. 
The K/-K misidentification probabilities arc fixed from 
the samples of the favored modes. 

The BB background populates the entire AE region. 
This background is fitted with an exponential PDF that 
models the tails of the backgrounds from the modes 
B~ -> D*n-, B~ ->■ Dp-, and B~ -)• D*K~ , which 
peak in the negative AE region, as well as combinatorial 
backgrounds. The NB PDF is obtained from BB MC 
samples, in which all known B and B meson decays are 
included. The charge asymmetry for the BB background 
is floated in the fit. 

The qq background distribution in AE is modeled by a 
linear function. The PDF of NB is obtained from a side- 
band sample of data: 5.20 GeV/c 2 < M hc < 5.26 GeV/c 2 
and 0.15 GeV < AE < 0.30 GeV. The charge asymme- 
try for the qq background is fixed to zero. 



The results of the fits of the above PDFs to the final 
data samples are shown in Fig. [T] The signal yields and 
the reconstruction efficiencies are listed in Table Q] Note 
that the rare charmless b —> s decay B~ — s- K + K~tt~ 
can peak inside the signal region for B~ — > [K + tt~]dK~ 
and be included in the signal yield. To estimate its contri- 
bution as well as contributions from B~ — > [K + K~]dtt~ 
and B~ -+ [K'n + ] D K~, we fit the M(Ktt) data side- 
bands: 1.815 GeV/c 2 < M(Ktt) < 1.845 GeV/c 2 and 
1.885 GeV/c 2 < M{Kir) < 2.005 GeV/c 2 . The side- 
bands arc chosen to avoid the contribution from B~ — > 
[K + K~]e>it~ caused by K/tt misidentification. We ap- 
ply the same fitting method used in the signal extraction 
to the sideband sample to obtain an expected yield of 
— 1.9lJ;5 even ts. For B~ — > [K + tt~]dtt~ , we also apply 
the requirement M(Ktt) < 1.915 GeV/c 2 for the side- 
band sample to avoid B~ — > [tt + tt~]eitt~ background, 
and obtain — 3.2 jlg'^. We do not subtract these back- 
grounds from the signal yields but instead include the 
errors on the yields in the systematic uncertainties. 
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TABLE I: Signal yields, reconstruction efficiencies and signif- 
icances. Charge-conjugate modes are included. The uncer- 
tainties shown are statistical only. 



Mode 



B~ 
B~ 
B~ 
B~ 



[K+tv-] d K- 
[K~n + } D K- 
■ [K + tt-] d tv- 



Yield 

T57T 
14.2 



Efficiency (%) Significance 



56.0 
3394 



-1+68 
-69 

165.0t&i 

1+245 
-244 



49164Z 



33.6 ± 0.4 
33.2 ± 0.4 
36.5 ±0.4 

35.7 ±0.4 



4.1a 



9.2a 



where the contributions to the systematic uncertainties 
are listed in Table [TTJ The uncertainties due to the 
AE PDFs for the DK signal, the Dtt signal, and the 
Dtt feed-across are evaluated by varying the shape pa- 
rameters by ±1(7. Those due to the DK feed-across 
are obtained by varying the width and the mean by 
±10%, which is the difference observed between the data 
and MC samples for the Dtt feed-across. The uncer- 
tainties from the NB PDFs for the DK and Dtt sig- 
nals (the Dtt feed-across) are estimated by obtaining 
PDFs from the region 0.01 GeV < \AE\ < 0.02 GeV 
(0.01 GeV < \AE - 0.05 GcV| < 0.02 GeV). Those due 
to the DK feed-across and the BB background are esti- 
mated by using the DK signal PDF. Those due to the qq 
background are estimated by using (Mbc AE) different 
sidebands. The uncertainties due to the K/tt misidenti- 
fication probabilities for the feed-across backgrounds are 
obtained by varying their values by their ±lcr errors. The 
uncertainty due to the charge asymmetry of the qq back- 
ground is obtained by varying it by ±0.02 (±0.005) for 
DK (Dtt), which is the uncertainty in the favored DK 
(Dtt) signal. A possible fit bias is checked by generat- 
ing 10,000 pseudo-experiments. The uncertainties in de- 
tection efficiencies mainly arise from MC statistics and 
the uncertainties in the particle identification efficien- 
cies. The total systematic uncertainty is calculated by 
summing the above uncertainties in quadrature. 



FIG. 1: AE (NB > 0.9) and NB (\AE\ < 0.03 GeV) dis- 
tributions for [K + tt~]dK~ (left) and [K + tt~]dtt~ (right). 
Charge-conjugate decays are included. In these plots, 
[K + tt~]dK~ components are shown by thicker dashed curves 
(red), and [K + tt^]dtt^ components are shown by thinner 
dashed curves (magenta). BB backgrounds are shown by 
dash-dotted curves (green) while qq backgrounds are shown 
by dotted curves (blue). The sums of all components are 
shown by solid curves (black). 



From the signal yields in Table |H we obtain 
K DK = [1.63i^(stat)t°;°3(syst)] x 10~ 2 



[3.28±° : 3 3 «(stat; 



-0.13V 
+0.12 
-0.18 



(syst)] x 10~ 



(5) 
(G) 



TABLE II: Summary of the systematic uncertainties. We use 
the notation "■ • •" to denote negligible contributions. 



Source 



Rl 



Al 



Al 



PDFs of AE 

PDFs of NB 

K/tt misidentification 

Asymmetry of qq background 

Fit bias 

Peaking backgrounds 
Efficiency 

Detector asymmetry 



~F2TT 
-1 



% 
-3.0 /0 

±0.2% 

+0.80/ 
-0.9 /o 

-1.1% 

-6.6% 
±1.7% 



+ 1.3c/ 
-1.2 /o 

±3.1% 
±0.0% 
±0.1% 
-0.5% 
-4.2% 
±1.5% 



Total 



+4.4o/ 
-7.8 /o 



TT7 



% 



±0.01 ±0.00 
+0.02 +f) nl 
-o.oi ^u-u 1 

±0.00 ±0.00 
±0.01 ±0.00 
-0.01 -0.00 
±0.03 +0.00 

±0.02 ±0.00 

+0.04 



hi). 02 
-0.01 



The significances oITZdk and TZd-k are estimated using 
the fit likelihoods by convolving asymmetric Gaussians 
denoting the systematic uncertainties Q, and listed in 
Table HI The significance for TZdk is 4.1cr, which consti- 
tutes the first evidence for the suppressed DK decay. 

The AE projections are shown separately for each 
charge of the B candidate in Fig. [2] We obtain 



Ad if 

•A-Dtt 



-0.391 



0.28 



(stat): 



-0.04 ± 0.11(stat)+^(syst), 



+0.02/ 



(7) 

(8) 



where the systematic uncertainties are evaluated in a sim- 
ilar manner as that done for TZdk and TZd-k (see Ta- 
ble UTJ) . The uncertainty due to the yield of the peaking 
backgrounds is obtained by varying the signal yield in the 
denominator of the CP asymmetry. The uncertainty due 
to the asymmetry of the peaking backgrounds is negligi- 
ble [2(|. To account for possible bias due to the charge 
asymmetry of the detector, we take the uncertainty in the 
asymmetry of the favored signal as a conservative limit 
on this effect. 
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FIG. 2: AE distributions (NB > 0.9) for [K+tt^dK- (left 
upper), [K~tt + ]dK + (right upper), [K + tt~]dtt~ (left lower), 
and \K~ ■k + ]dk + (right lower). The curves show the same 
components as in Fig. [T] 



Assuming Eqs. ©-(gj) and r B = 0.1 [15J, the val- 
ues of TZdk and Adk are restricted to the ranges [0.2, 
2.5] X 10~ 2 and [—0.9, 0.9], respectively, depending on the 
values of </>3, <5g, and 5d- Our results are consistent with 
these expectations. The small experimental uncertainties 
in TZdk and Adk thus provide important additional in- 
formation on (j)^. The experimental results for TZd-k and 
Adtt are also consistent with the standard model [151 . 



In summary, we report measurements of the suppressed 
decay B~ — > [K + ir~]Dh~ (h = K,ir) using the full 
772 x 10 6 BB pair data sample collected with the Belle 
detector. We use a neural network-based method [lj| to 
discriminate qq background from signal, impose a D* ± 
veto, and employ a two-dimensional fit to extract the sig- 
nal. These steps, along with a 20% increase in the data 
sample, result in a significant improvement compared to 
the previous analysis ||. We obtain the first evidence 
for a DK signal with a significance of 4.1<r, and report 
the most precise measurements to date of the CP asym- 
metries and ratios of the suppressed decay rate to the 
favored decay rate. Our results will provide important 
ingredients in a model- independent extraction of </>3 . 
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